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Abstract 

We give a complete description of the genus expansion of the one-cut solution to 
the generalized Penner model. The solution is presented in a form which allows us in 
a very straightforward manner to localize critical points and to investigate the scaling 
behaviour of the model in the vicinity of these points. We carry out an analysis of 
the critical behaviour to all genera addressing all types of multi-critical points. In 
certain regions of the coupling constant space the model must be defined via analytical 
continuation. We show in detail how this works for the Penner model. Using analytical 
continuation it is possible to reach the fermionic 1-matrix model. We show that the 
critical points of the fermionic 1-matrix model can be indexed by an integer, ra, as it 
was the case for the ordinary hermitian 1-matrix model. Furthermore the ra'th multi- 
critical fermionic model has to all genera the same value of ~/ str as the ra'th multi-critical 
hermitian model. However, the coefficients of the topological expansion need not be 
the same in the two cases. We show explicitly how it is possible with a fermionic matrix 
model to reach a ra = 2 multi-critical point for which the topological expansion has 
alternating signs, but otherwise coincides with the usual Painleve expansion. 
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1 Introduction 



The hermitian 1-matrix model with a polynomial potential is by now completely un- 
derstood. All critical points have been localized and classified. The different types of 
critical behaviour are indexed by an integer, m, and the m'th multi-critical model is 
characterized by ~f str taking the value — ^ [1]. Furthermore it is well known that if 
one considers a potential consisting of a linear plus a logarithmic term one can reach 
a critical point for which ^y str = and for which logarithmic scaling violations occur 
at genus zero [2]. This is the critical point of the Penner model [3]. The possibility of 
multi-critical behaviour for generalized Penner models, i.e. models with a logarithmic 
term and a higher degree polynomial was pointed out in references [4, 5, 6]. However, 
only the genus zero and the genus one behaviour was addressed. 

In reference [7] a complete description of the genus expanded 1-cut solution to 
the hermitian 1-matrix model with a generic polynomial potential was given. Here we 
generalize this description to the case where in addition a non polynomial term, namely 
a logarithm, appears in the interaction. The solution is presented in a form which 
allows us in a straightforward manner to localize all critical points and to investigate 
the scaling behaviour of the model in the vicinity of these points. All types of multi- 
critical points can be addressed and the analysis can be carried out for any genus. The 

outcome of our analyses is that no other values of ^f s tr than ^f s tr — or ^f s tr — ? 

m = 2, 3, . . . are possible. 

In some regions of its coupling constant space the generalized Penner model can 
only be defined via analytical continuation. We show in detail how this procedure works 
for the Penner model itself extending the analyses of references [4, 5]. By analytical 
continuation it is possible to reach the fermionic 1-matrix model [8]. As an application 
of our study of the generalized Penner model we show that the possible types of multi- 
critical behaviour for the fermionic 1-matrix model can, as in the hermitian case, be 
indexed by an integer, m. Futhermore the m'th multi-critical fermionic model has to 
all genera the same value of ~/ str as the m'th multi-critical hermitian matrix model. 
However, the coefficients of the terms in the topological expansion will not necessarily 
be the same for a m'th multi-critical point obtained from a hermitian and a fermionic 
matrix model respectively. For example it is possible in the fermionic case to find a 
m = 2 multi-critical point for which the topological series has alternating signs but 
otherwise coincides with the usual Painleve expansion. Having alternating signs, this 
series might very well be Borel summable. 
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2 The model and its loop equations 



The generalized Penner model is defined by the partition function 

Z = e N2F =[ #exp(-iVTr [U(<f>)]) (2.1) 

JNxN 

where the integration is over hermitian N X N matrices and 

00 a- 

U(<f>) = V(<f>)+tlog<f>, V^) = J2-^- (2-2) 

3=1 J 

We define expectation values in the usual way and introduce the generating functional 
for 1-loop averages by 

W'-(p)=<^Tr-^). (2 .3) 

With the normalization chosen above the genus expansion of the free energy and the 
1-loop correlator read 

00-1 00-1 

F = ^W^ W ^ = Y. Wg W g { P )- (2-4) 

3=0 3=0 

The loop equations of the model express the invariance of its partition function under 
field redefinitions. To derive the loop equations in an appropriate form it is convenient 
to consider the following transformation of the field, <f> 

^^ + e E^i = ^ + e (-r^\) ■ (2 - 5) 

„m? + \P{P-<P)J 

The reason why we do not include a n = term in the sum above is that we wish to 
avoid the appearance of terms of the type Tr </> _1 from the variation of log <f> in the 
action. Introducing the shift of </>, (2.5) in (2.1) we get to the first order in e 

T r (^ff - NT, (W,^) } = 0. ,2.6, 

We now introduce as usual corresponding to the matrix <f> an eigenvalue density p(X). 
With our interaction it is necessary to require that the support of the eigenvalue density 
does not cross the branch cut of log A. Then the equation (2.6) can be written as 

- f ^^^M + -w( P ) = (w( P )f + ^wn w ^ ^ 

p Jc Ziri p — u p l\ 2 aV(p) 

where d y^ is the loop insertion operator 



— = _vJ-A (2S] 

dV(p) ^ p^ dg 3 K ' ] 
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and where C is a curve which encloses the support of the eigenvalue density but which 
does not cross the branch cut of the logarithm, neither encloses the point u = p. It is 
easy to see that (2.7) can be written in the form 

c 2ttz p-u K ' v N 2 dV{p) KFJ K ' 

provided W(p) fulfills the following equation 

. ^-.V'(uj)W(uj) - tW(0) = 0. (2.10) 

For t = the validity of the relation (2.10) is ensured by the invariance of the partition 
function under field redefinitions of the type <f> — > <f> + e. For t ^ we can also 
formally derive the relation (2.10) with VK^(O) = (^Tr </> -1 } by considering the same 
transformation of the field. However, the ill defined quantity Tr </> _1 appears. By 
equating the ^ terms in (2.9) we see that a W(p) which fulfills (2.9) will automatically 
fulfill the requirement (2.10). Hence a solution of (2.9) is automatically a solution 
of (2.7). The loop equation (2.9) is of exactly the same form as the one of the hermitian 
1-matrix model without a logarithmic interaction term. In reference [7] an iterative 
procedure for solving the latter loop equation was developed and explicit results for 
W g (p) as well as F g for g = 1 and g = 2 were given. Furthermore the general structure 
of F g and W g (p) for any g was described. From the results of reference [7] we can read 
off the complete perturbative one-cut solution of the model (2.1). 

3 The genus zero solution 
3.1 General results 

With the assumption that the singularities of W(p) consist only of one square root 
branch cut (which in the case of W(p) being a solution of (2.9) is equivalent to p(X) 
having support only on one arc in the complex plane [9]) and with the normalization 
W(p) — > 1/p as p — > oo, the genus zero contribution to W(p) for the model (2.1) can 
be written as 

ww= i/iM v ' M+ -U''-iM,r (3-D 

VfV 2Jc2tti p-u \(u - x)(u - y) j K ' 

where x and y are determined by the boundary conditions 

, , f du {V'(lo) + ^} , s 

- Zc2^;_; )( „ _,,="• < 3 - 2 > 

B 2 (x,y) = i / V ; J = = 2 . (3.3) 

Jc 2ttz ^J( u _ _ y ) 
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It is easy to verify by direct calculation using (3.1) that the condition (2.10) is indeed 
satisfied for the genus zero solution due to the boundary conditions (3.2) and (3.3). 
With V(cf>) given by (2.2) the boundary equations can be written as a set of algebraic 
equations. For deg(V((f>)) = D the degree of the boundary equations is D if t = 
and 2D if t ^ 0. Hence for t = an explicit solution can be found for deg(V((f>)) < 4 
whereas for t / an explicit solution can be found only in the case of a linear or a 
quadratic potential. It is easy to see that for t = l,x = y = zis& possible solution of 
the boundary equations for any potential. Namely setting x = y = z in the boundary 
equations (3.2) and (3.3) we find 

B 1 (z,z)=V'(z)-- B 2 (z,z) = zV'(z)+t. (3.4) 

Hence with t = 1 and B\(z,z) = the boundary equation B 2 (z } z) = 2 is automatically 
fulfilled. Conversely it follows that the only way in which we can have x = y = z is by 
t being equal to 1. As noted in reference [4] for t = 1 we have Wo(p) = ^ independently 
of V(cf>). This result follows easily from formula (3.1). Another case for which some 
information can be extracted from the boundary equations even for a generic potential 
is the case t = 2. For t = 2 one finds that if the potential V{cf>) is of a definite parity x 
and y will lie on the imaginary axis. If V{cf>) is even x and y will both be situated either 
on the negative or on the positive imaginary axis. If V{cf>) is odd one has x = —y. By 
solving the boundary equations we determine the branch points of W(p). However, 
this does not fix the position of the cut of W(p) } i.e. the support of the eigenvalue 
distribution. The missing information is encoded in the function G(X) [10] 

G(X) = I" (U'(p)-2W (p))dp. (3.5) 
Jy 

The support of the eigenvalue distribution is an arc connecting y to x along which 
G(X) is purely imaginary and which is embedded in a region where Re(Cr) < 0. We 
are used to considering a model of the type (2.1) as meaningful if the branch points x 
and y are real and the cut lies along the real axis. There are regions of the coupling 
constant space where such a situation can not be realized. In some of these regions it 
is still possible to attribute a meaning to the model by analytical continuation. The 
eigenvalue density is defined according to the rules given above and is now a curve 
in the complex plane. We note that one might obtain by this prescription a complex 
valued partition function. A lot of information about the model (2.1) can be extracted 
without knowing the precise location of the branch points or the cut of W(p). For 
instance all critical points can be easily localized and classified. All we have to assume 
is that the branch cut of W(p) does not cross the branch cut of the logarithm. When 
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working with contour integrals this assumption turns into the assumption that the 
contour does not cross the branch cut of the logarithm. 

By means of the solution (3.1) we can calculate the genus zero contribution to the 
susceptibility, Xo 

* = -i R (:5 - 6 » 

By differentiating (2.1) one gets 

X = dt { N Tr l0g = to 2™ l0g " ~dT ■ (3 - ?) 

From the expression (3.1) one obtains by direct differentiation using Eqs. (3.2) and 
(3.3) 

/ \ 1/2 

dW (p) f did 1 / (to — x)(lo — y) \ 



dt Jc 47rz (p — io) to \ (p — x)(p — y) J 

= 1 f 1 l i Xy 1 ) (38) 

2\p p\(p-x)(p-y)J (ip-x)(p-y)) 1/2 f 

Substituting into Eq. (3.7) and compressing the contour C to the cut of the logarithm, 
we get 

_i r d [i U *y V /2 i 1 

X ° 2io P \p p \(p - x)(p - y) ) ((p_ x )(p_ y ))i/2j 

1 f (2(xyy/ 2 - (x + y)Y 

We stress that this formula holds for a generic potential. 

3.2 The Penner model 

Here we show how the procedure of analytical continuation works in a simple case, 
namely the case of the Penner model 

g x = -1, g t = 0, i > 1. (3.10) 

For this model the boundary equations read 

B 1 (x,y) = -l-—L- = 0, (3.11) 
\xyyi 1 

B 2 (x,y) = -±( x + y )+t = 2 (3.12) 
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and we note that we recover from (3.9) immediately the well known result [4] 



Xo = -log(l -t)- log* 



(3.13) 



Solving (3.11) and (3.12) for x and y gives 



x = (t - 2) + 2(1 - t) 1/2 , 

y = ( t _ 2 )-2(l-t) 1 / 2 . 



(3.14) 
(3.15) 



The analyticity structure of W(p) and hence the location of the support of the eigen- 
value distribution is quite different for different values of t. A detailed analysis of the 
analyticity structure based on the function G(A), which can be calculated explicitly, is 
carried out in Appendix A. Here we just summarize the results. 

t < 0: The branch points x and y are real and negative and the cut of W(p) is the 
straight line connecting x and y. This position of the cut of W(p) implies that 
(xy) 1 ! 2 = l^j/l 1 / 2 which is in accordance with equation (3.11). The branch cut of 
the logarithm can for instance be placed along the negative real axis (figure 1). In 
the limit t — > the cut of W(p) collides with the cut of the logarithm (figure 2). 

< t < 1: The branch points x and y are again situated on the negative real axis. 
However, the cut of W(p) now looks as in figure 3. We note that the presence 
of the loop implies that (xy) 1 / 2 = — Ixy] 1 / 2 which is again in accordance with 
equation (3.11). For t — > 1 the branch points x and y collide (figure 4). 

t > 1: The branch points x and y get an imaginary part and are each others complex 
conjugate. The cut of W(p) is an arc which connects y and x and which intersects 
the real axis on its positive part (figure 5). In particular for t = 2 the model (2.1) 
can be interpreted as a fermionic matrix model [8]. We will return to this point 
in section 6. 

4 The solution to higher genera 

Explicit solutions for W\(p) and W2(p) as well as the general structure of W g (p) for 
g > 1 can be read off from reference [7] and it is easy to show that W(p) fulfills 
the requirement (2.10) to all genera. From the analysis of reference [7] we also have 
a detailed knowledge of the free energy of the model. To express the higher genera 
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Figure 1: The support of the spectral density for t < (the bold line) and the branch 
cut of the logarithm (the thin line). 
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Figure 2: The support of the spectral density for t 
cut of the logarithm (the thin line). 



(the bold line) and the branch 
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Figure 3: The support of the spectral density for < t < 1 (the bold line) and the 
branch cut of the logarithm (the thin line). 
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Figure 4: The support of the spectral density for ^ — > 1 where x — > y (the bold line) 
and the branch cut of the logarithm (the thin line). 
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Figure 5: The support of the spectral density for t > 1 (the bold line) and the branch 
cut of the logarithm (the thin line). 



contributions to F (as well as to any multi-loop correlator) it is convenient to introduce 
instead of the coupling constants ({</;}, t) a set of moments {Mfc, J^} 

Jc 2vrz {to - x) 1 / 2 {to - j,)*+i/2 - v > 

The genus one contribution to F reads 

Fx = -— In Mi - — In J x - -In J. (4.3) 

24 24 6 y ' 



For g > 1, F g takes the following form. 

F g = J2 {ai... a s ;P 1 ...p l \a,l3,j) g M ajP^ ' 9 ~ ' ^ ^ 



where d = x — y } the brackets denote rational numbers and a, f3 and 7 are non-negative 
integers. The indices «i, . . . , a s , . . . , fii take values in the interval [2, 3g — 2] and the 
summation is over sets of indices obeying the following restrictions. 

s-a<0, ( s -a = 0)^(s = a = 0) (4.5) 



9 



l-/3<0, (I - (3 = 0) ^ (I = (3 = 0), (4.6) 
(a-s) + (/3-l) = 2g-2, (4.7) 

- 1) + - 1) + 7 = % - 4, (4.8) 

i=l j=l 

5-l<7<4^-4. (4.9) 

We note that contour integrals like (4.1) and (4.2) are calculated by taking residuals 
at zero and at infinity. Hence in our formalism we effectively do not need to know the 
precise location of the support of the eigenvalues, although it is of course essential in 
order to get a complete understanding of the model. We see that for real x and y } F g 
will automatically be real while in the case of models with x and y complex, which 
must of course be defined via analytical continuation, F g will in general be complex. 
In the following we investigate the coupling constant space of the one-cut solution to 
the model (2.1) defined when necessary and possible via analytical continuation. We 
localize and classify its critical points and show how in the vicinity of any of these one 
can define a continuum theory using a double scaling prescription. In particular we 
show that there exist critical points for which x and y are complex but for which we 
obtain real valued F g } s for all g. 

5 Classification of critical points 

This section is devoted to the study of the critical properties of the one cut solution 
to the model (2.1). Actually the most straightforward way to localize the critical 
points is by taking a glance at the expressions (4.3) and (4.4) for the higher genera 
contributions to F. It is obvious that there are several circumstances under which 
the higher genera contributions to the free energy become singular. One possibility 
is that Mi or J\ or possibly both acquire a zero of some order. This was the only 
possibility in the case of the usual 1-matrix model and corresponds to the situation 
where a number of extra zeros accumulates at one or possibly both ends of the support 
of the eigenvalue distribution. However with the logarithmic term in the action we 
also have the possibility of d becoming zero. As shown in section 3 the only way in 
which d can become zero is by t being equal to one. Of course we can also have a 
situation where both d = and and Mi (=Ji) acquire a zero of some order. A third 
possibility for singular behaviour is that some of the higher moments diverge as we 
approach a certain point in the coupling constant space. This will be the case if at 
the critical point x — > (or y — > 0). For x — > we must require that t — > in order 
that E>i(x } y) is well defined. However, the critical point t = can be approached in 
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a number of different ways which are all characterized by x = 0{t). That t = and 
t = 1 are critical points of the model (2.1) was first noted in reference [4] where the 
genus zero contribution to the susceptibility was analyzed. Here we will carry out a 
complete perturbative analysis of the scaling behaviour of the model in the vicinity of 
its critical points addressing all types of multi- critical behaviour. We will approach the 
critical points by fixing the coupling constants at their critical values, {gf}, and 
scaling only t. As appears from the discussion above the critical points of the model 
are naturally divided into three types. 

Type 1: t$ {0, 1}, M x = M 2 = . . . = M m _ x = 0, M m ^ 0, ra > 2, 
Ji = J i = ■ ■ ■ = J n -i = 0, J n ^ 0, 1 < n < m. 

Type 2: / 1. \/, \l . . . . M m _ x = 0, M m ^ 0, m > 

Type 3: / 0. .W„ ,\/, ... M m _ x = 0, M m ^ 0, m > 0, J = J x = . . . = 

J n - X = 0, J n ^ 0, 1 < n < m 

Here and in the following it is understood that for moments of the type Mk({gi}, t = 0) 
and Jk({gi},t = 0), t should be set to zero before the contour integration is carried 
out. We note that since we have understood that all scaling is associated with t the 
usual ra'th multi-critical points of the hermitian 1-matrix model do not immediately 
fit into this classification. For the usual ra'th multi-critical points one uses a scaling 
prescription where t is kept fixed at t c = and the coupling constants {g{\ are scaled. 
As it will appear in the following for t {0, 1} both scaling prescriptions lead to the 
same critical behaviour. For t c = 1, keeping t fixed at t c does not make sense (cf. 
to (4.4)). To characterize the different types of critical behaviour let us denote by A# 
the renormalized version of t (i.e. A# ~ it — t c )) and let us introduce the critical index 
Istr by 

Fg „ A g- 7str )(i- 3 ) 5 g > h (51) 
5.1 Type 1 critical points 

Keeping the coupling constants fixed at their critical values and assuming for 

a given value of t the branch points of W(p) to be x and y we find expanding the 
boundary equations B x (x } y) and B 2 (x } y) keeping only leading order terms 

c m (x-x c ) m M c m + c n (y-y c ) n J c n = (5.2) 
c m {x-x c ) m x c M c m + c n {y-y c ) n y c J c n = (t c - t) (5.3) 
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where 

(2m — 1 V 

=M m ({tf},f c ), J c k = J k ({ 9 n,t c ), c m = { J- (5.4) 

ml Z m 

If we introduce a renormalized coupling constant A R by 

(t c -t) = a m A R (5.5) 

we see that 

(x - x c ) ~ aA R /m , (y - y c ) ~ {a m A R ) 1 l n (5.6) 

while of course d ~ a since per assumption t ^ 1. Furthermore expanding the moments 
one finds 

M k ~ (aA R lm ) m -\ ke[l,m], J/ ~ (a m A#)^, le[l,n]- (5-7) 

(We note that in order for (5.7) to be true it is necessary that neither x nor y tend 
to zero at the critical point. However, as pointed out above x — > (or ?/ — > 0) is not 
possible for t finite.) From the expression (3.9) for Xo it follows that the leading non 
analytical term in F is given by 

2+ — 

F = const ■ A R m (5.8) 

where const is a non universal constant. (We remind the reader that we have assumed 
that m > n.) Furthermore from the relation (5.7) it follows that in the scaling limit 
the genus 1 contribution to the free energy develops a logarithmic singularity of the 
following type 

1 f in — 1 n — 11 

Fi = < 1 >\ogA R 1 < m, 1 < n < m. (5.9) 

24 I m n J 

A given term of F g} g > 1 will scale with a negative power of a, P£, given by 

s / 
777 777 

Pi = a( m -l)+f3(m--)-Y(m-a i )-Y(m-f3 J -) (5.10) 

= (g - l)(2m + 2) - 7 + |(/3 - /) + - 1) J (1 - ^) (5.11) 

where to obtain the second equality sign we have made use of the relation (4.8). Now 
it follows from (4.5), (4.6) and (4.9) that for n < m the dominant terms of F g are those 
for which 

13 = 1 = 0, i = g-l (5.12) 
These terms hence do not depend on any J-moments and have 

s 

( a _ s ) = 2 5 -2, £(<*,-- 1) = 3s - 3 (5.13) 



i=i 
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In deriving (5.11) we implicitly assumed scaling for all M moments involved. By 
a slight modification of the argument given above it is easy to convince oneself that 
terms which contain a moment Mk with k > m will be subdominant in the limit a — > 0. 
Hence in this limit F g can be written as 

P V- / I \ M ai ...M as 

F g = ^ ( ai ...a s \a) g (5.14) 

l<aj<m IV1 1 Uc 

All terms in (5.14) scale as given by 

Max{P 9 _} = (g - l)(2m + 1) (5.15) 

and we see the possibility of a double scaling limit emerging. Furthermore, bearing in 
mind (5.5), we recover the behaviour (5.1) with 

Istr = -1/m. (5.16) 

For n = m all terms with 7 = g — 1 contribute in the scaling limit and F g becomes 
a sum of two terms of the type (5.14), one which involves Af-moments and one which 
involves J- moments (cf. to reference [7]). Obviously the relations (5.15) and (5.16) 
hold also in this case. The type 1 multi-critical points of the generalized Penner model 
clearly belong to the same universality class as the usual ra'th multi-critical points of 
the hermitian 1-matrix model. 

5.2 Type 2 critical points 

Let us consider a critical point ({ <?«•}, t) = ({gf}, 1) for which Mi = M 2 = . . . = M m _i = 
0, M m ^ 0, m > 1 and let us approach the point by setting ({</;}, t) = ({g 8 c },t). We 
note that Mj: = Jj: since x c = y c . Assuming for a given value of t the branch points of 
Wo(p) to be x and y we find by expanding E>i(x } y) keeping only leading order terms 

m 

J2 M c m Ci c m _ t {x - x c y (y - y c ) m ~ l = (5.17) 

8 = 

which tells us that (x — x c ) ~ (y — y c ). The expansion of B 2 (x } y) leads to the following 
relation 

m + l 

J2 {*c M c m+1 + M m ) Ci c m+1 _ 8 (x - x c y (y - y c ) m+1 ~ l =t c -t (5.18) 

8 = 

(The vanishing of the terms of order (x — x c ) m is due to (5.17).) We now define a 
renormalized coupling constant A# by 

t c -t = a m+1 A R (5.19) 
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and get 

d ~ ak]l {m+1) (5.20) 
The moments scale in the following way 

M k ~ J k ~ (aA)/ (m+1) ) m " fc , A; G [0,m]. (5.21) 



Using the expression (3.9) it is easy to see that the singular behaviour at genus zero is 
caused by the scaling of d and that 

F = ^-A^logA^ (5.22) 

m + 1 

as was also found in reference [4]. Furthermore using (4.3), (5.20) and (5.21) we see 
that in the scaling limit we have independently of m 

F 1 = -^\ogK R . (5.23) 

This property of the multi-critical Penner models was conjectured (but not proven) in 
reference [6]. Here we can address also the question of the behaviour of these models 
for g > 1. A given term of F g} g > 1 will scale with a negative power of a, P 9 _ given by 

P£ = (( a _ s ) + (^_/))( m _i) + ^( a _i) + ^_i) +7 (5.24) 

i=i j=i 

= (flf-l)(2m + 2) (5.25) 

Here we find the amazing result that all terms in our expression (4.4) are potentially 
relevant for the scaling limit. By potentially relevant we mean relevant for m sufficiently 
large: We have in deriving (5.24) and (5.25) implicitly assumed that all moments scale. 
A closer analysis reveals that as in the previous case terms which contain a moment M^, 
k > m can be neglected in the limit a — > 0. Now bearing in mind the relation (5.19) 
we see that independently of m we have 

Istr = 0. (5.26) 

Among the type 2 critical points we find the critical point of the Penner model, g\ = — 1, 
g\ = 0, i > 1. 1 It is well known [2] that this model exhibits the scaling behaviour 
characteristic of a theory describing 2D quantum gravity interacting with matter with 
central charge, c = 1. However, it appears that this scaling behaviour occurs for all 

1 For the Penner model one actually has to keep next to leading order terms in the expansion of 
the boundary equations. However, it its easy to show that the conclusions stated above hold also in 
this case. 
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type 2 critical points with M£ ^ 0. Let us note that lor these models all scaling of the 
higher genera contributions to the free energy is due to the scaling of d and the terms 
which are important in the limit a — > are those for which 

7 = 4(^-1), s = l = 0, a + /3 = 2(g-l) (5.27) 

In particular these terms do not depend on any moments other than Mi and J\. 

5.3 Type 3 critical points 

For simplicity let us consider a point for which J\ ^ 0. The generalization to the case 
J\ = J2 = • • • = Jn-i = 0, J n ^0,l<n<mis straightforward although more tedious 
than in the case of the type 1 multi critical models. (We remind the reader that for 
moments of the type Mk({gi},t = 0), Jk({gi},t = 0) it is understood that t should be 
set to zero before the contour integration is carried out.) Again we keep the coupling 
constants {g{\ fixed at their critical values and assume at a given value of t the branch 
points of W(p) to be x and y. After performing a few rearrangements we can write the 
boundary equations (3.2) and (3.3) as 

c m x m (xy c ) 1 t 2 MZ l + c 1 (y-y c )JZy c = (5.28) 
c m x m {xy c ) x l 2 M c m + t = (5.29) 

Hence if we introduce a renormalized coupling constant A# by 

t = a m+1 / 2 A R (5.30) 

we have 

x~aA™, (y-y c )~a m+1 / 2 A R (5.31) 
Making use of the expression (3.9) for Xo we nn d 

^0 = ^— A^logA^. (5.32) 
Am + 1 

For the moments we have 

M fc ~a m " fc , k £ [0, 00] (5.33) 

while the J-moments do not scale. Hence from (4.3) we see that the genus 1 contribu- 
tion to the free energy behaves in the scaling limit as 

F 1 = -—( m ~\ ) logA^. (5.34) 
24 Vm + 1/2; 5 y J 
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To handle the higher genera contributions to F we note that except for the fact that 
in the present case all M- moments scale the behaviour (5.33) is the same as that of 
a type 1 multi-critical model with (m,n) = (ra, 1). Hence we have the same value of 
max{Pf!_ } as in that case which, bearing in mind the relation (5.31) tells us that 

Istr = 0. (5.35) 



6 The fermionic 1-matrix model 

In this section we consider the following 1-matrix model 

Z = e N2p = J «$exp (-iVTr (6.1) 

where and ^ are N X N matrices whose matrix elements are independent Grassmann 
variables. We take the interaction to be of the following type 

oo ( . 

= ( 6 - 2 ) 

k=i k 

For this model we define the 1-loop correlator by 

^H^Tr^) (6.3) 

and the genus expansion of the free energy and the 1-loop correlator looks as in equa- 
tion (2.4). To derive the loop equations of the model it is convenient to consider a 
transformation of the fields (vp, *&) of the following type 

* ->■ * + e— ^— , (6.4) 

# -> (6.5) 



Under such a transformation the measure changes as 



\p-\jf\ff I I \P(P- **) 



(6.7) 



2 In case we were considering in stead of the matrices \?) with Grassmannian matrix elements 
complex matrices (<f>\<f>) a similar shift of variables would lead to the following change of the mea- 
sure [11] 

^-^(l + ^Tr (—L^))^. (6.6) 

What causes the difference between this formula and formula (6.7) is the fact that in comparison with 
the complex matrix model the fermionic matrix model has an additional factor (—1) appearing in the 
change of the measure for transformations of the type \? — ► \? + eM^^M*)™ , n > 1. 
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Using the result (6.7) it is easy to show that the invariance of the partition function 
under the transformation (6.4), (6.5) leads to the following loop equation 

1 1 dU uV ' {u) -W(u) + 2 W( P ) = (W( P ) f + ^^W(p) (6.8) 



p Jc 2iri p — to p N 2 dV(p) 

where the notation is as in section 2. This equation is identical to equation (2.6) for 
t = 2 [8]. Hence we learn that when analytically continued to t = 2 the hermitian 
matrix model (2.1) is equivalent to the fermionic matrix model (6.1). In particular we 
know from section 3 and section 4 the complete perturbative one-cut solution to the 
model (6.1). 

The boundary equations for t = 2 reveal an interesting feature of the fermionic 
model. Inserting t = 2 in (3.2) and (3.3) we find 

We see that the roles of the boundary equations are in a sense opposite to what was 
the case for the ordinary hermitian 1-matrix model. The concepts symmetrical and 
non symmetrical are interchanged, i.e. an odd potential leads to the branch points of 
W(p) being placed symmetrically with respect to the origin while an even potential 
leads to a non symmetrical placing of the branch points. Let us consider the case of 
an odd potential, i.e. 

oo f . 

n**) = E^zr(**) 2 * +1 - I 6 - 11 ) 

In this case, if we set x = —y = z the boundary equation (6.10) is trivially fulfilled, 

since there is no residue at infinity. The remaining boundary equation reads 

oo 2 

E^ic J ^- ( -^ = (6.12) 

which admits a solution of the type z = ifi. In order to recover in the limit <?2j+i — > 0, 
j > the previously obtained results for the Penner model we should take the negative 
square root of (— z 2 ) 1 ! 2 . However, if we want to recover in the limit <?2j+i — > 0, j > 
a gaussian model we should take the positive square of (— z 2 ) 1 ! 2 . In the following we 
will consider the latter possibility since we wish to have an interpretation of our model 
in terms of random surfaces. In any case we see that for a fermionic matrix model 
with odd potential we have for the branch points of W(p) } x = —y = ifi. This has the 
following implications for the moments 

M n = (-l) n J n . (6.13) 
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It is obvious that the fermionic model possesses only critical points of type 1. To 
investigate the scaling behaviour in the vicinity of these points we will set = {g-gf} 
and scale g to 1. As mentioned earlier for type 1 critical points the critical behaviour 
is the same whether the scaling is associated with t or g. Hence the critical points of 
the fermionic model belong to the same universality class as the critical points of the 
usual hermitian 1-matrix model. The value of ~/ str for the ra'th multi-critical fermionic 
model coincides with the value of ~/ str for the ra'th multi-critical hermitian model to 
all genera. The critical indices being identical does not mean that the topological 
expansion in the vicinity of a ra'th multi-critical point is the same for the fermionic 
and the hermitian 1-matrix model, however. For the fermionic 1-matrix model the 
branch points of W(p) } i.e. x and y will in general be complex and hence F g will in 
general be a complex quantity. In the case where the potential of the fermionic model 
has a definite parity it is easy to investigate how the branch points being complex 
affects the topological expansion. 

Let us to begin with consider the case of an odd potential. Inserting x = —y = ifi 
in the definitions (4.1) and (4.2) we find 

M n = -M n , .1 = ' ./. (6.14) 

where M is real and given by 

" J 2iri (lo — fj) n+1 l 2 (io + fj) 1 ! 2 (-//)"// 1 ' 

and J n = ( — l) n M n (cf. to equation (6.13)). In addition we of course have 

d = 2ifi = id (6.16) 

where d is real. Inserting the relations (6.14) and (6.16) into the expression (4.4) we 
find for a given term f g of our expression (4.4) for the genus g contribution to the free 
energy 

f g = ( z )«+^-e; =1 °>i-T,U p—<j g ( 6 .i7) 

where f g is real and appears from f g in (4.4) by replacing M , J and d by M , J and d. 
Using the relations (4.7) and (4.8) one finds that 

F g = (-ly-'F,. (6.18) 

Now it seems as if a possibility of a topological expansion with alternating signs 
emerges. To test the viability of this possibility we must of course study in detail 
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the Af's. Let us do so restricting ourselves to the vicinity ol a m = 2 multi-critical 
point. 

At this point it is instructive to recall what the situation was like for the hermitian 
1-matrix model without the logarithmic term. To be specific let us consider the model 
where only g\ and g 3 are different from zero. For this model only Mi and M 2 (and J\ 
and J 2 ) are non vanishing (where we now understand that the moments are defined 
without the t term). The support of the eigenvalue distribution (in the perturbative 
region of the coupling constant space) is asymmetric and we denote it as The 
model has a m = 2 multi-critical point, namely 

9i = 93 = -I- (6-19) 

Let us approach this critical point by setting gi = g ■ gf. Expanding the boundary 
equations in powers of (x — x c ) and (y — y c ) we recover the equations (5.2) and (5.3) 
with (m,n) = (2, 1). In particular inserting (5.2) in (5.3) we find 

(1 - g) = c 2 (x - x c fM c 2 . (6.20) 

In the scaling limit all dependence on J-moments disappear and for the Af-moments 
we have 

M 2 = M 2 C = ^<0, (6.21) 
M x = ij .U, (x - x c ) > 0. (6.22) 

Furthermore for d we have 

d = x c — y c = d c > 0. (6.23) 

Now let us return to the fermionic model with odd potential and let us try to 
localize a m = 2 critical point, i.e. a point for which M = Mi = 0. If we assume for 
simplicity that only g\ and g 3 are different from zero we get the following equations for 
9i, 9z and // c 

M C = -I^ + 5l > c -2 = 0, (6.24) 
Ml = ^-2 = 0. (6.25) 

It is easy to see that the following solution is possible 

93 = h 5i c =2^, /^c = 2 1 / 3 . (6.26) 

As before we approach the critical point by setting gi = g ■ gf and expand our (single) 
boundary equation in powers of z — z c (where z = ifi). After a few rewritings where only 
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leading order terms are kept we can represent our boundary equation in the following 
form 

(g-l) = c 2 M c 2 n c {n-n c f (6.27) 

which looks completely similar to the corresponding equation of the ordinary hermitian 
1-matrix model (cf. to (6.20)). In this case we find in the scaling limit for the moments 
of interest 

M 2 = -G/ 3 c + 4) = -253<0, (6.28) 
M x = M 2 C (// - fi c ) = -2g c 3 (/i - fi c ) > 0. (6.29) 

Furthermore we of course have 

i=2/i c >0. (6.30) 

Since Ml, M 2 and d have the same signs as Mi, M 2 and d for the hermitian model 
we get for the fermionic m = 2 multi-critical model a topological expansion with 
alternating signs. Apart from the signs the topological expansion of the fermionic 
model is identical to the usual Painleve expansion (for suitable normalization of the 
cosmological constant). The question of the possible change of signs in the topological 
expansion of the higher multi-critical fermionic models requires further investigation. 
However, on the basis of the above presented example we put forward the conjecture 
that for all ra'th multi-critical points of the fermionic 1-matrix model with odd potential 
we have in the scaling limit for the genus g contribution to the free energy, F^ erm 

pferm _ ^_-^y-l pherm 

where F^ erm is the genus g contribution to the free energy of a ra'th multi-critical 
hermitian model obtained from a symmetrical potential. 

Let us close this section by commenting on the fermionic matrix model with an even 
potential. In this case an analysis of the boundary equations shows that the branch 
points of W(p) are again imaginary. However, now they both lie on either the positive 
or the negative real axis. (If the potential does not have a definite parity the branch 
points will in general have a real as well as an imaginary part.) Also for the symmetric 
fermionic matrix model we get a set of relations like (6.14). However, the prefactor ^ 
is replaced by ^Vr- This modification results in (6.18) being replaced by 

F g = F g (6.32) 

Analyzing in detail the ra = 2 multi-critical point is in this case a lot more complicated 
than in the case of an odd potential, because the cut is no longer symmetric and the 
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moments involve more terms due to the larger degree of the potential. Of course the 
analysis of the higher ra'th multi-critical points is even more difficult. However, we put 
forward the conjecture that the topological expansion in the vicinity of a ra'th multi- 
critical point obtained from a fermionic matrix model with even potential coincides 
with the topological expansion in the vicinity of a ra'th multi-critical point obtained 
from a non symmetrical hermitian matrix model. 

7 Conclusion 

The combination of the loop equations and the moment description constitute a pow- 
erful tool for studying higher genera contributions in matrix models. In particular the 
moment variables make it very simple to localize critical points and to investigate the 
scaling behaviour in the vicinity of these points. 

Previously the moment technique was developed only for 1-matrix models with 
polynomial potentials. Here we have extended the method to the case where in addition 
a non polynomial term, namely a logarithm, appears in the interaction, i.e. to the 
generalized Penner model. The investigation of the coupling constant space of the 
generalized Penner model revealed no unknown types of scaling behaviour. Our solution 
of the generalized Penner model provided us with a solution of the fermionic 1-matrix 
model and the possible types of critical behaviour of this model turned out to be 
characterized by the same value of ^y str as the ra'th multi-critical points of the usual 
hermitian 1-matrix model. A feature of the fermionic 1-matrix model which deserves 
some attention, though, is that it has a ra = 2 multi-critical point for which the 
topological expansion has alternating signs but otherwise coincides with the usual 
Painleve expansion. Having alternating signs this series might very well be Borel 
summable. It would be interesting if the fermionic nature of the model could be given 
a world sheet interpretation. 

One interesting prospect of the present work is the possibility of generalizing the 
moment technique to the hermitian two-matrix model which is known to be capable 
of describing all minimal conformal models coupled to gravity [12]. The study of 
loop equations for the hermitian two-matrix model has been initiated [13], but explicit 
results for correlators are few and limited to genus zero. Another promising application 
of our results is the genus zero analysis of matrix models in dimensions, D > 1 [14]. 
Here the task would be to localize critical points and devise possible continuum limits. 
The genus zero loop equations were obtained in reference [15]. Of particular interest 
is the hermitian matrix model with a logarithmic potential for which these equations 
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can be exactly solved for any D [16]. 



Appendix A The eigenvalue support for various t 

Let us investigate in some detail the question of the location of the support of the 
eigenvalue distribution for the Penner model 

U(<f>) = -<f> + t\og<f>. (A.l) 

From equation (3.14) and (3.15) we know explicitly the position of the endpoints x 
and y of the distribution as a function of t. As explained earlier the exact location of 
the support can be found only by studying the function G(X) defined by (3.5). The 
support is an arc connecting y to x along which G(X) is purely imaginary and which 
is embedded in a region where Re(G(A)) < 0. Using our solution (3.1) we can write 
down an explicit expression for G(X). We find 

G(X) = - [ x l(p- x y/*(p- y y/* 

h P 

= ~ [(P ~ x) 1/2 (p ~ yf' 2 + (2 - t) log (2 \(p - xfl 2 (p - yf' 2 + p + (2 - f )] ) 

2l-t(p- x y/i(p-yy/i + p(2-t) + t*] Y\ x 
p J 

/ J y 

where we have on the way made use of the boundary equations (3.11) and (3.12). 

Below we have shown the sign variation of Re(G(A)) in some spectacular cases 
assuming the cut of the square root (p — x) 1 / 2 ^ — j/) 1 / 2 to be the straight line connecting 
x and y. Before proceeding to analyzing these plots let us note certain characteristica 
of the function Re(G(A)) which can easily be read off from (A. 2). First we see that as 
p — > ±oo we have Re(G(A)) — > =poo and as p — > we have Re(G(A)) — > +oo (provided 
t ^ 0). Futhermore it is easy to convince oneself that if x and y are real Re(G(A)) = 
for A £ [j/,^]. However, if x and y are imaginary we have Re(G(A)) ^ for A's on the 
imaginary axis between x and y. 

Figure 6 shows the case t = 1. We see that the branch cut [y,x] has collapsed to 
a point and that there is a loop encircling the point A = along which Re(G(A)) = 0. 
Inside the loop Re(G(A)) > 0. The appearance of such a loop going through the point 
x is a common feature for all t < 1, t ^ 0. For t < the support of the eigenvalue 
distribution is simply [j/,x]. However, for < t < 1 with the cut of W(p) only reaching 
from y to x we have actually not fulfilled the requirement that we should take the 
negative square root of (xy) 1 / 2 (cf. to section 3). However, this requirement will be 
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fulfilled if we deform our original cut so that it includes also the loop. Then we will 
have Re(G(A)) < both inside and outside the loop and hence the loop belonging to 
the support of the eigenvalue distribution. When t > 1 the points x and y become 
complex and each others complex conjugate. In figure 7 we have plotted Re(G(A)) for 
t = |. We see that there is still an arc connecting x and y along which Re(G(A)) = 0. 
This arc passes the real axis to the right of the origin. (We note that the discontinuity 
of Re(G(A)) along the vertical line connecting x and y of course has its origin in the cut 
of the square root (p — xY^ 2 (p — y) 1 / 2 being placed along this line). As in the case of 
< t < 1 we can now deform our original cut into the arc with Re(G(A)) = connecting 
x and y, thereby identifying this arc with the support of the eigenvalue distribution 
and fulfilling the requirement concerning the sign of (xy) 1 / 2 . Finally in figure 8 we have 
plotted Re(G(A)) in the fermionic case t = 2. The figure examplifies the statement 
that the roles of even and odd potentials are interchanged for the fermionic 1-matrix 
model. We see for the odd potential V{cf>) = —cf> a completely symmetric picture. The 
analysis of the cut structure is left to the reader. 
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Figure 6: The sign variation of Re(G(A)) for t = 1. Re(G(A)) < in the gray region 
and Re(G(A)) > in the white region. 



25 




26 



-4 -2 2 4 



Figure 8: The sign variation of Re(G(A)) for t = 2. Re(G(A)) < in the gray region 
and Re(G(A)) > in the white region. 
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